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The Quenched Approximation in Health and in Sickness * 

Claude Bernard'^and Maarten Golterman 

'^Department of Physics, Washington University 
St. Louis, MO 63130-4899, USA 

We present results for physical quantities computed in quenched chiral perturbation theory and compare them 
with the corresponding unquenched expressions. We also point out an apparent theoretical problem of the 
quenched approximation. 



1. Introduction 

We have recently proposed a systematic la- 
grangian approach to chiral perturbation theory 
(ChPT) for the quenched approximation of QCD 
[1,2]. This approach is based on the observa- 
tion that QCD can be quenched by adding a 
number of scalar (ghost) quarks to the QCD la- 
grangian with exactly the same quantum numbers 
and masses as the physical, fermionic quarks, so 
that the combined determinants cancel exactly 
[3] . This defines a lagrangian formalism for QCD 
in the quenched approximation. The extended 
lagrangian possesses a graded [/(3|3)l(XiC/(3|3)_r 
chiral symmetry. 

One can then develop ChPT for this larger 
symmetry group in the usual way, assuming the 
existence of a Goldstone meson multiplet which 
corresponds to the breakdown of this group to 
its diagonal subgroup. This multiplet does not 
only contain the usual Goldstone mesons, but 
also fermionic mesons which are bound states of 
a physical quark and a ghost anti-quark or of 
a physical anti-quark and a ghost quark, and 
ghost mesons which are ghost quark anti-quark 
bound states. One has to take into account the 
graded nature of the symmetry group, by replac- 
ing traces and determinants by super-traces and 
super-determinants in the construction of the ap- 
propriate invariants. 

In fact, the symmetry gets broken down 
to the semi-direct product group [SU{3\3)l <Xi 
SU{S\3)r](s)U{1) by the anomaly,^ and corre- 



spondingly there exists a "super- ry'", $0, which 
is invariant under this reduced symmetry group. 
$0 is a linear combination of the r/' and its ghost 
partner, fj': 

$0 = str$ = -^(77' - f]'), 

where str denotes the super-trace and $ is a 
graded hermitian 6x6 matrix describing the 
Goldstone meson multiplet (for details, see ref. 
[2]). 

All these considerations lead to an 0{p^) chiral 
lagrangian (in a notation similar to the notation 
used in ref. [4]) 

L = ^Voi^o) + 14($o)str(9^Sa^S^) 

where the functions Vi can be chosen to be real 
and even, S = exp(2i$//), and 
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where m„, and are the quark masses. 

The special role of the 77' in the quenched ap- 
proximation becomes clear from considering the 
quadratic part of L for the 77' and 77' fields which 
(for degenerate quark masses) reads 
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direct product in refs. [1,2], 
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This leads to a propagator in the rj-f]' sector 
1 / 1 



mt 



+ 



1 / (p^ - rri^y 



1 1 
1 1 



,(1) 



where fi^ = ^Vq'{0) is the parameter which in the 
full theory gives the singlet part of the rj' mass. 
(We have set V5 = for simplicity. One can show 
that vertices coming from V5 contribute only at 
higher order in a combined expansion in 1 /Nc and 
M to quantities considered in this talk.) 

Because /x^ appears in the numerator of eq. 
(1), it is clear that the 77' docs not decouple for 
large /x^ in quenched ChPT, unlike the case of 
full QCD. This can also intuitively be understood 
from the "quark flow" approach. For a discussion 
and applications of the quark flow approach, see 
refs. [5-7,2]. 

2. Comparison between quenched and full 
ChPT 

In full ChPT the ratio fK/fn at one loop for 
m„ = rud = m is [4] 

loop 
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where - j^^n.p,^^ ^2 
coefficient of an 0{p'^) term. In quenched ChPT 
the result is [1,2] 

1— loop 
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167r2/2 
+ {ms - m)L5^q^. 

These two expressions can however not be com- 
pared directly, since in general the parameters in 
the full and quenched ChPT lagrangians will not 
be equal: Ls^qu ^ is.fuU- This implies that in or- 
der to compare full and quenched ChPT we will 
have to consider physical quantities independent 

of bare parameters. In the full theory, such a 
1 _i 

quantity is frj-n fj^^^ [4]. However, in quenched 
ChPT with nis ^ m, the inherits the r]' double 
pole (eq. (1)) through mixing. 



A similar physical quantity can be defined in 
Nf = A QCD with quark masses = md = m 
and m.s = nis' = m' . The quantity we wish to 
consider here is 



R- 
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where tt, K and tt' denote the du, su and s's 
mesons, respectively. This quantity is natural to 
compute in a quenched simulation, since one has 
the quark propagators for various masses. The 
total number of flavors, Nf, is irrelevant, as only 
the valence quarks play a role in the calculation 
of any quenched quantity. In full ChPT, 
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whereas in quenched ChPT, 
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where we used the tree level relation m|- = 
i(m^ + m^,). Note that the logs in the full 
and quenched expression have completely differ- 
ent forms, stemming from their different physical 
origins. 

To get a feeling for the numerical difference we 
have substituted two sets of numerical values for 
the quantities appearing on the right-hand side of 
eqs. (2,3). First, we have considered the "real" 
world, with m„ = 140 MeV and = 494 MeV 
(with TOtt' computed from its tree level relation 
to m^r and mx), and with fjl^/Z = (500 MeV)^ 
from the full QCD 77' mass. Secondly, we have 
used the quark and meson masses from a typical 
simulation (ref. [8], /3 = 6.2), in particular m = 
0.007, m' = 0.03 and a"^ = 2.9 GeV. (Again we 
took ^2/3 = (500 MeV)2.) For the "real" world 
we get 



iJfuii = 1.023, R, 



■qu 



1.066 



(a difference of 4.2%), whereas for the /? = 6.2 

data we find 

i?fuu = 1.022, i?qu = 1.014 
(a diflFerence of 0.8%). 
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We have also looked at an Nf = 3 quantity 
independent of the bare parameters (using results 
from [2]), namely 
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{uu) m]^o - 1^1+ {uu) ' 



(In QCD, the divergences of {qq) cancel in x [4]-) 
At tree level, the values for full and quenched 
ChPT coincide, but again, at one loop 

Ax = x'-'°°P - X*'"''^ 
differs. 

For the real world, we find 

^tree ^ q gg2 

and 

Axfuu = -0.002, Axqu = -0.064 

(the difference between Xfuii and Xqu is 6.3%). 
Prom the /? = 6.2 data of ref. [8] with quark 
masses m„ = 0.01, rrid = 0.02 and rrig = 0.03 we 
obtain 



Xtree — 



1 

2' 



Axfuii = -0.017, Axqu = -0.026 

(a difference of 1.9% in x)- 

The one-loop corrections for the full and 
quenched quantities are unrelated for both R and 
X (for an example where they are the same, see 
ref. [6]). However, for both examples (R and x), 
the difference between the full and quenched case 
is small. This is due simply to the fact that the 
one-loop corrections themselves are small. Note 
that the differences between the quenched and 
full case grow with decreasing quark mass, due to 
the fact that the quenched one-loop corrections 
grow with decreasing quark mass. We will return 
to this point in the next section. To our knowl- 
edge, the results presented in this section consti- 
tute the first controlled calculation of the differ- 
ence between quantities in the quenched and full 
theories. 

We have been using an estimate for /z^ from 
the full theory, whereas of course its quenched 
value need not be equal to the full value. It would 



therefore be nice if we could measure /Xq„ from 
quenched data. One could in principle determine 
/iqy from a fit of the pion mass to the form [2] 
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Wc have tried to fit both the /? = 6.2 and /? = 6.0 
data from ref. [8] (not only for TOtt but also for 
{'tptp)). However, we have not been able to un- 
cover any chiral logs in these data; one apparently 
needs smaller errors and lower quark masses. 

3. Problems with quenched QCD? 

Let us again consider the quenched result for 
fK/f-TT (cf. sect. 2). Rewritten in terms of quark 
masses the result is 
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It is clear that the limit my, with fixed, 
which is well defined in the full theory, does not 
exist here [1,2]. Stated in other words, taking 
both TO„, nis 0, one is left with a result de- 
pendent on the ratio rris/mu. (A similar prob- 
lem in (V'V') has been discussed in ref. [6]. See 
also ref. [7].) This seems to imply that no chi- 
ral limit exists for quenched QCD! We note here 
that the bare quark mass appearing in the chi- 
ral lagrangian is analytic in the bare quark mass 
which appears in the QCD lagrangian (in both 
the full and the quenched cases). This means 
that redefinitions of the quark mass [6] cannot 
solve the problem. 

The problem arises from IR divergences due to 
the double pole term in the rj' propagator (cf. eq. 
(1)). The first question one may therefore ask 
is whether this double pole term is really there, 
or whether it might be softened by summing the 
quenched perturbative expansion. 

This can be phrased as follows (for simplicity 
we consider the case of degenerate quark masses) . 
In the summed theory, the parameter /Xq^, appear- 
ing in eq. (1) will become momentum dependent. 
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(In the quenched theory we may take eq. (1) 
as a definition of /iqu(p) to all orders. Similarly, 
in the full theory l^f^nip) is defined to all orders 
as the singlet part of the r]' self-energy.) If now 
Hquip) ^ for p — > 0, the double pole would be 
softened, leading to a less divergent chiral limit. 
It appears however, that this will not occur: 
We first observe that to leading order in 1/Nc 

and secondly, that 
mLi(p = 0)^0 

in particular when the quark mass m is set to 
zero. The latter statement follows from the fact 
that the rj' is a well-behaved meson in full QCD 
with a nonvanishing mass in the chiral limit. As- 
suming the validity of the 1/Nc expansion, we 
conclude that 

4uiP = 0) 7^ 0, 

which means that the double pole is a true feature 
of the quenched theory! 

This argument does not prove that no resum- 
mation could ameliorate the problem. On this 
issue, however, we have the following remarks: 

• Sharpe has summed a class of diagrams in the 

case of degenerate quark masses [6,7]; this 
did not lead to a less singular result. 

• In the nondegenerate case there are many more 

relevant diagrams, and it is not clear that a 
systematic resummation can be carried out. 

• Any resummation method would have to work 

for each quantity that does not have a chi- 
ral limit at one loop. Many such quanti- 
ties exist. If some resummation works for 
each quantity, one would expect that there 
is a general explanation of why it works. 
We have shown above that this explanation 
cannot be a softening of the rj' double pole 
term, which nevertheless seems to be the 
origin of the problem. 

In view of these remarks, we believe that the 
nonexistence of a chiral limit is a real problem of 
the quenched theory, and not an artifact of chiral 
perturbation theory. 



4. Conclusions 

We have calculated in a controlled way the dif- 
ference between the quenched and full versions of 
QCD for two physical quantities. Such calcula- 
tions can be extended straightforwardly to many 
other physical quantities. 

We have also pointed out that quenched QCD 
does not seem to have a well-behaved chiral limit, 
due to the peculiar role of the r]' in the quenched 
approximation. We believe that this is a real 
problem of the quenched theory, and not an arti- 
fact of ChPT. In order for the quenched approx- 
imation to remain a key tool in lattice QCD, it 
will be essential to gain a better understanding of 
this problem. 
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